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m : Abstract 

o 
o 

We discuss a ID many-body model of distinguishable particles with local, momen- 
Cp , tum dependent two-body interactions. We show that the restriction of this model 

^ I to fermions corresponds to the non-relativistic limit of the massive Thirring model. 

' This fermion model can be solved exactly by a mapping to the ID boson gas with in- 

verse coupling constant. We provide evidence that this mapping is the non-relativistic 
limit of the duality between the massive Thirring model and the quantum sine-Gordon 
model. We also investigate the question if the generalization of this model to distin- 
guishable particles is exactly solvable by the coordinate Bethe ansatz and find that 
this is not the case. 



1. Introduction. In this paper we present, discuss, and solve a non-relativistic many-body 
5^ \ system of particles moving in one space dimension (ID) and interacting with a particular 
local, momentum dependent two-body potential. As we will explain, this model is the 
natural fermion- analog of the ID boson gas. 

The ID boson gas is one of the famous exactly solvable many-body models. It describes 
non-relativistic bosons moving in ID and interacting with delta-function two-body interac- 
tions, and it was solved by Lieb and Liniger a long time ago [1] (a nice textbook discussion 
of this model and its solution can be found in Chapter I of reference [2]). We mention 
only in passing the considerable recent interest by experimental physicists triggered by a 
proposal of an experimental realization of this model in [4]. 

As mentioned, the particles in the ID boson gas interact via a delta-function interaction. 
Due to the Pauli principle, this kind of interaction is trivial for fermions, and thus interesting 

^Corrected version of an article with the same title in J. Phys. A: Math. Gen. 37 No 16 (2004), 4579- 
4592. 
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fermion models with such an interaction require additional internal degrees of freedom [5, 6]. 
Our ID many-body model is without internal degrees of freedom and with a particular 
local, translation invariant interaction which is non-trivial for fermions. It is defined by the 
following Hamiltonian {d^j = d/dxj), 

N 

j=l j<k 

with an arbitrary number N of particles moving on the real line, —oo<Xj<oo (we will also 
mention some generalizations of our results to an interval of length L with periodic or anti- 
periodic boundary conditions, < Xj < L); the real parameter A determines the coupling 
strength. Note that the interactions depend not only on the particle distance Xj — Xk but 
also the momentum difference pj — = —1(9^;. — d^i^)- As we will see (Paragraph 2), 
due to this the Pauli principle is circumvented: this interaction is non-trivial on fermion 
wave functions, while it is trivial on boson wave functions. We will derive this model as 
non-relativistic limit of the massive Thirring model [7] (Paragraph 3), in the same way 
as the boson gas can be obtained as non-relativistic limit of 0^-theory in 1+1 dimensions 
(see Appendix B.2). We find that this fermion model can be solved exactly by mapping 
it to the ID boson gas with the coupling replaced by its inverse (Paragraph 4; as we will 
explain, this result is equivalent to the duality observed previously in [8]). This relation 
between our fermion model and the ID boson gas is reminiscent to the famous duality 
between the massive Thirring model and the quantum sine-Gordon model [9], and we will 
present arguments that it actually is the non-relativistic hmit of the latter (Paragraph 5). 
A natural question is if the generalization of this model to distinguishable particles remains 
exactly solvable by the coordinate Bethe Ansatz, but, as we shortly discuss in Paragraph 6, 
this is not the case. 

Since the massive Thirring model is known to be integrable (in certain formal meanings 
of this word), it is perhaps not too surprising that its non-relativistic limit in equation (1) is 
exactly solvable. It thus is worth recalling that, despite of various interesting partial results 
[10, 11], the Thirring model has not been solved in full detail. It thus is interesting that its 
non-relativistic limit can be solved and studied by the much simpler methods which have 
been developed for the ID boson gas. 

In our derivations of non-relativistic limits in Paragraph 3 we start with the formal 
definition of the quantum massive Thirring model, perform expansions in 1/ (mass x c) with 
c the velocity of light, and we use physical arguments to justify our ignoring certain terms. 
In this way we arrive at a non-relativistic model which is well-defined, in the same spirit as 
reference [12]. It should be possible to make this procedure mathematically precise using 
the method proposed in [13]. As mentioned, the arguments in Paragraph 5 are somewhat 
heuristic. The other results are mathematically precise. We tried to keep the main text 
short, but for the convenience of the reader we included two appendices: In Appendix A 
we give a complimentary treatment of the singular interaction in our model. Appendix B 
contains details of our non-relativistic limits. 

Remcirk: In the first version of this paper ( which unfortunately has already appeared in J. 
Phys. A: Math. Gen.) we presented an argument which seems to show that the generalized 
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model defined by he Hamiltonian in equation (1) and for distinguishable particles is exactly 
solvable by the coordinate Bethe Ansatz [3]. This argument was incorrect. 

The mistake leading to our wrong conclusions was a subtle point in Yang's arguments 
[3], and we were mislead by one of our sources to use an unfortunate notation covering up 
this point. To be more specific: while equation (C4) in the published version of this paper 
is not incorrect, it is misleading the way it is written, and our interpretation spelled out 
in Appendix C.l, Remark 2, led us astray: the twist we proposed to save the validity of 
the Yang-Baxter relations leads to an inconsistency at another point which we missed. The 
correct interpretation of this equation is now given in section 6. 

It is worth mentioning that we discovered this mistake when we tried to generalize our 
result to other types of singular interactions, checking the validity of the Bethe Ansatz for 
the three particle case directly using symbolic computer programs written in MAPLE and 
MATHEMATICA (to be sure two of us wrote independent programs) . While these programs 
nicely confirmed Yang 's result for the delta interaction case, they showed that the coordinate 
Bethe Ansatz for our model is consistent if and only if the eigenfunction either has boson 
or fermion statistics. 



2. Two particle case. To get a physical understanding of our model it is instructive to 
first consider the two-particle case N = 2. Introducing x = xi — X2 and ignoring the trivial 
center-of-mass motion, H in equation (1) reduces to the following simple Hamiltonian, 

h^-dl + AXdJ{x)d,, (2) 

whose eigenfunctions x{^): 2; e R, are defined by satisfying {dl + E)x{x) = for x 7^ and 
the following boundary conditions, 

x'(o+) - x'(-o+) = 

x(0+)-x(-0+) = 4Ax'(0+), (3) 

with the prime indicating differentiation. Indeed, these are the boundary conditions ob- 
tained by integrating hx = Ex twice: first from x = —0"*" to x > where x'(0) interpreted 
as the average of the left- and right derivative, and then once more from x — —0^ to 0"*" (it 
is instructive to verify this formal argument by checking that the solutions below indeed 
satisfy hx — Ex)- The solutions of this are obtained by simple computations, 

X+{x) = cos{kx) 

. . sm(kx) / X /, X /.X 

X-W = 2AA; ^ ^g^(^) cos(A;a:) (4) 

with corresponding eigenvalue E — k"^. For real k these all are scattering states, and 
for A < there is one additional bound state for k — i/2A with energy E — — 1/4A^. 
Thus positive and negative values of A correspond to the repulsive and attractive cases, 
respectively. As already mentioned, the boson wave function x+ is unchanged by the 
interaction, while the fermion wave function X- is modified, opposite to what happens for 
the delta-function interaction. It is worth noting that, in converting the interaction in 
equation (2) into the boundary conditions in equation (3), we have used a regularization 
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procedure which consistently avoids divergences which would occur in a naive treatment of 
this singular interaction (this is explained in more detail in Appendix A). 

In a similar manner one finds that the eigenfunctions x of the Hamiltonian in equation 
(1) for arbitrary N are given by the solutions of (Yl - d^^ +E)x{xi, . . . , xn) = in all regions 
of non-coinciding points, together with the following boundary conditions 

(^Xj 9xi^x\xj=Xk+0+ i^^^j ^Xk^X\xj=Xk—0'^ 

X\xj=Xk+0+ ~ X\xj=Xk-0+ — '^^(^(^Xj ~ dx^^x\xj=Xk-0+ (5) 

(we used that d^^-xk — i'^xj — dx^)/'2). It is straightforward to check that these bound- 
ary conditions are trivially fulfilled for all non-interacting boson eigenfunctions x+ — 
S^PeSjv ^-^P^Sj i^Pj^j)- They are, however, non-trivial for fermions. 



3. Non-relativistic limit of the massive Thirring model. We now derive the non- 
relativistic limit of the massive Thirring model [7] and show that it is identical with the 
second quantization of the many-body Hamiltonians in equation (1). The Thirring model 
can be (formally) defined by the quantum field theory Hamiltonian Ji — Tio + TYint where 
the free part is the usual Dirac Hamiltonian in ID, 

«.^/d.:(.Ul)(-'t7^° ,.a:- E„) it) (^) 

with m > the fermion mass, and the interaction is 

Hint ^^gjdx: V'iV+V'-V'- : (7) 

(see e.g. equation (2.1) in [10]) with g the coupling constant and the dots indicating normal 
ordering; the ip^^^ = 'iIj^^\x) are fermion field operators obeying the usual canonical anticom- 
mutation relations (CAR) {i^^ix) , ip^^{y)} = 6{x — y) etc., and Eq is a parameter allowing 
us to change the reference energy which we will fix later to a convenient value. One can 
diagonalize T^o by Fourier transformation and diagonalization of a 2 x 2 matrix, which corre- 
sponds to a particular canonical transformation (■0+,'0-) ~^ (^+) (see Appendix B.l). 
We expand in powers of 1/mc and obtain, in position space, 

^, = -L(*,±*_^ + (8) 

and Tio = ^0 + with Hq = ± j dx : ^^[(mc^ ^ Eo) - dl/2m + . . :, where the 
dots are for higher order terms in 1/mc. The positive- and negative states of the non- 
interacting model are now decoupled, and it is straightforward to compute the interaction 
in terms of the new fields '^±. To obtain the non-relativistic limit we set Eq = itk? and 
assume that mc^ is large. In this case we can ignore the negative energy degrees of freedom 
the non- interacting ground state is such that all the negative energy states are filled 
and the positive energy states empty (Dirac sea), and due to the large energy gap 2mc^ the 
interactions involving the filled states, in particular those across the gap, can be neglected if 
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one is only interested in the low-energy physics. We thus drop all terms in the Hamiltonian 
involving the fields and in leading non-trivial order in 1/mc we obtain the following 
Hamiltonian, 

^non-rel ^ Jdx (-dl)^ + : ((^.^t) (a^^)^t^ _ *t (fi,^^) (gl^^t) ^) . (9) 

with ^ = ^+ obeying CAR and annihilating the non-interacting vacuum, ^|0) = 0; we 
used : [^'^(a;)\E'(x)]^ := 0, i.e., the lowest order term vanishes due to the Pauli principle, and 
thus the leading non-trivial interaction involves derivatives. It is straightforward to verify 
that this non-relativistic quantum field Hamiltonian TYnon-rei is the second quantization 
of our many-body Hamiltonian H in equation (1): for 2m = 1 and g/{2mc)'^ — —A, the 
eigenvalue equation T^non-reil-^) — E\N) for A^-particle states 

\N) = j d''xx{x,,...,XN)^\x,)---^\xNm (10) 

is equivalent to Hx = Ex- Note that A < corresponds to > 0, in agreement with what 
one should have expected from the fact that the massive Thirring model has bound states 
for > (see Eq. (2.15b) jf in [10]), whereas the sign of A is such that the attractive case 
corresponds to A < (see Paragraph 2 above). 



4. Solution I: Fermion model. We now determine all fermion eigenfunctions x of 
the Hamiltonian in equation (1). Due to the fermion statistics we only need to determine 
X — xixi,X2, ■ ■ ■ , xn) in the fundamental wedge 

Aj : xi < X2 < . . . < xpf. (11) 

For the same reason, the boundary conditions in the first line of Eq. (5) are automatically 
fulfilled, and the ones in the second line simplify to 2x\xj=Xk+o+ = 2A((9a;^. — dxf.)x\xj=x^+o+ 
where we only need to consider the cases j — k + 1. Thus the equations determining our 
eigenfunctions are d^. + E)x — and 

(^x.+i - dx^ - ^)xU,+i=x^+o+ = 0. (12) 

Comparing with equations (2.1a), (2.4a) in [1] we see that these conditions are identical with 
the ones determining the eigenfunctions of the ID boson gas defined by the Hamiltonian 

N 

Hb^-Y. 9% + E '^(^^- - (13) 

3=1 j<k 

at coupling 

CB^J (14) 

in the fundamental wedge A/. Since the latter eigenfunctions are well-known, we can im- 
mediately write down all eigenfunctions of our model 

x{xi,X2, ...,xn)= TT ( A[a^ - dx^] + 1 ) det [exp(i/cja;fc)] (15) 

\ / l<i,k<N 

l<k<j<N ~ ~ 
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in A/, and the corresponding eigenvalues are E — A;| (this exphcit formula is apparently 
due to Gaudin [14]; see Chapter I in [2]). 

In this paper we restrict ourselves to particles moving on the full line, but it is interesting 
to note that many of our results can be extended to the finite interval of length L, < 
Xj < L, with periodic or anti-periodic boundary conditions, 

X(xi, . . . , xjv) = e'^xi^i + L,...,xn): (16) 

and similarly for all other arguments Xj, with 77 = or tt. Similarly as for the ID boson 
gas this yields the following conditions for the allowed momentum values, 

(these are the so-called Bethe equations; see e.g. Chapter I in [2]). Comparing with the 
Bethe equations for the ID boson gas (equation (2.2) in [2]) we sec that the duality above 
remains true for finite interval if we choose in our model periodic boundary condition 
[f] = 0) if is even and anti-periodic boundary conditions (77 = vr) if is odd. In 
the thermodynamic limit L, N ^ 00 such that p = N/L remains finite the difference 
in boundary conditions becomes irrelevant, and thus all thermodynamic properties of our 
model are the same as the known thermodynamic properties of the ID boson gas [16] at 
inverse coupling, = 1/A. It would be interesting to know if there are any observables 
which can distinguish these two models. 



5. Non-relativistic limit of the quantum sine-Gordon model. We now present 
evidence that the relation of our fermion model to the ID boson gas found above is the 
non-relativistic limit of the duality between the massive Thirring model and the quantum 
sine-Gordon (qSG) model [9]. In the main text we will argue that the qSG model reduces 

to (f)\^i-ih.eoiY for large (effective) mass. The result then follows since 0f^^-theory in the 
non-relativistic limit is identical with the second quantization of the ID boson gas (the 
details of this latter part of the argument are deferred to Appendix B.2). 

The qSG model can be formally defined by the Hamiltonian TisG = 'Hq + Tif with the 
usual free boson Hamiltonian 

n^ = \jdx: (c^n^ + ct>[-dl + {mcf\ct>) : (18) 

and the interaction 

= /dx :^[l-cos/?0]-^0^: (19) 

with boson fields = (j){x) = (p^ and their conjugate variables 11 = dt(f)/c^ obeying the 
usual canonical commutation relations (CCR), [n(x), (f){y)] — —i5{x — y) etc.; a and (3 are 
coupling parameters. It is important to note that, while the bosons in the qSG model are 
massless, the interaction generates a mass m with 

{mcf = a. (20) 
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We moved this mass term to the free part of the Hamiltonian so that the Taylor series of 
the interaction starts with the forth order term, 



In the non-relativistic hmit we get, in leading order 1/mc, 

1 ($ + $t + . . .) (22) 

where are non-relativistic boson fields obeying the CCR [$(a;), $^(|/)] = S{x — y) (see 
Appendix B.2). Thus the coefficient in front of the n-th order term in the interaction is 
oc m^~"/9^"~^c^, suggesting that, if the mass is large, one only needs to take into account 
the leading term n = 2 of the interaction. We thus conclude that, for large values of a, the 
qSG model has the same non-relativistic limit as 0^_,_^-theory. Using that we find that the 
qSG Hamiltonian, in leading orders of 1/m and 1/c, reduces to 

<on-rei = / ^^Kx){-dlmx) - ^ : ^\xMx)^\xMx) : (23) 

where normal ordering is defined with respect to the vacuum |0) obeying $|0) = (sec 
Appendix B.2 for more details). This Hamiltonian now is well-defined, and for 2m = 1 and 
(/3c/4)^ = —cb it is identical with the second quantization of the ID boson Hamiltonian in 
equation (13) (see equations (1.1)-(1.12) in [2]). 

In Paragraph 4 we found a duality between the fcrmion Hamiltonian defined in equation 
(1) and the boson Hamiltonian in equation (13) with the relation of coupling parameters 
given in equation (14). It is interesting to compare this to Coleman's duality between the 
massive Thirring model and the qSG model (see equation (1.9) in [9]), 

Inserting the relations A = ~-g/(? and cb = —{(3c/A)^ which wc obtained in the non- 
relativistic limits in Paragraph 3 and above, we obtain the relation in equation (14) up to a 
factor 7r^/4 (the 1 on the r.h.s. in equation (24) disappears in the limit c — > oo). We regard 
this agreement up to a numerical factor of order one as strong evidence that the the duality 
found in Paragraph 4 is indeed the non-relativistic limit of Coleman's duality (note that an 
exact agreement cannot be expected since we ignore the rcnormalization of parameters in 
the qSG and massive Thirring models [9]). Note, however, that this argument only applies 
to the attractive case A < 0, whereas the duality in equation (14) is true also for A > 0. 

It is important to note that Coleman's duality provides also an identification of field 
operators in the qSG and the Thirring models (see equations (1-10) and (1-11) in [9]), but 
we do not see how this identification appears in our non-relativistic limits. We therefore 
regard the arguments in this paragraph only as a heuristic explanation of the duality in 
equation (14). It would be interesting to substantiate it in greater depth. 
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6. Generalization of the model to distinguishable particles. We now discuss the 
generalized model with the Hamiltonian in equation (1) but for distinguishable particles. 
We follow Yang [3] and make the following Bethe Ansatz for the eigenfunctions, 

X = ^ Ap(Q) exp (iE^Li kpjXQj) for xqi < xq2 < . . . xqn, (25) 
PeSN 

for all Q e S'jv, which implies E — It is straightforward to adapt Yang's computation 

to our boundary conditions in equation (5). Inserting the Bethe Ansatz we obtain 

i{kp, - kp^,+i))[ApTXQTi) - Ap{QTi)] = i{kp, - kp^,+i))[Ap{Q) - AptXQ)] (26) 
Ap{QT,) + AptXQT.) - Ap{Q) - AprM = '2XKkp, - kp^,+^))[Ap{Q) - AprM] (27) 

where Tj is the transposition interchanging i and i + 1. This implies 

(1 - iX[kp^i+i) - kpi])Ap{Q) = ApT.iQTi) - iA[A;p(,+i) - kpilAprM- (28) 

We now use that the permutation group Sn acts on the coefficients Ap{Q) by the regular 
representation it! — > i2 as follows, 

Ap{QR) = {RAp){Q) ^ Rq,Q'MQ') (29) 

with Rq,q' — Sq'^qr. We thus can insert ApTi{QTi) — (TiApTi){P) and write this latter 
relation as follows, 

Ap = Yi{kp^i+i) - kpi)ApTi (30) 
where Ap here is a vector with N\ elements Ap{Q) and 

As explained by Yang [3], consistency of the Bethe Ansatz is equivalent to Yi{u) satisfying 
certain consistency conditions. The first one (unitarity) is fulfilled in our case, 

Yi{-u)Yi{u) = f, (32) 

but the second one (Yang-Baxter relations) is not: 

Yi{v)Yi+^{u + v)Yi{u) ^ Yi+^{u)Yi{v + u)Yi+^{v). (33) 

We conclude that the generahzed model with distinguishable particles is not exactly solvable 
by the coordinate Bethe Ansatz. 



7. Final comments. It is well-known that, in addition to the delta-function interac- 
tion which has been studied extensively in the context of intcgrable many-body systems, 
there are other local interactions which are physically very different [17]. Recently it was 
found that one particular such interaction leads to an exactly solvable many-body system 
of fermions in ID which has a remarkable duality to the ID boson gas [8]. In this paper we 
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found a natural physical interpretation of this fermion model: we showed that the boundary 
conditions used to define the model in [8] naturally arise from the iV-body Hamiltonian in 
equation (1) which describes particles with local, momentum dependent two-body interac- 
tions. We also showed that this Hamiltonian arises as non-relativistic limit of the massive 
Thirring model, and we argued that the above-mentioned duality to the ID boson gas comes 
from the well-known duality of the Thirring model to the quantum sine-Gordon model. We 
then proposed a generalization of this model where the particles are distinguishable, but 
we found that, different from the delta interactions case [3], this model cannot be solved 
exactly by the coordinate Bethe Ansatz. 

As discussed in Chapter 1.4 of [17], quantum mechanical point interactions in ID leading 
to the boundary conditions in equation (3) have been studied extensively in the literature 
from a different point of view, and apparently it has been interpreted as a (^'-interaction 
(see [18]). The interpretation we give in this paper is very different and, as we hope to have 
convinced the reader, more natural. 

We believe that our results show that, from a physical and mathematical point of view, 
the model defined in equation (1) is equally interesting as the delta-function interaction 
model given in equation (13). It thus would be worthwhile to explore this model further, 
e.g., extend our results to the finite interval with suitable boundary conditions etc. 

Acknowledgments. E.L. would fike to thank Jens Hoppe and Antti Niemi for asking 
stimulating questions, Jouko Mickelsson for helpful comments, and Martin Hallnas for 
checking some of our formulas. E.L. would also like to thank the Erwin Schrodinger Institute 
in Vienna for hospitality where this work was started. CP. was supported by Deutsche 
Forschungsgemeinschaft under the Emmy-Noether programme. E.L. was supported in part 
by the Swedish Science Research Council (VR) and the Goran Gustafsson Foundation. 

Appendix A: Physical interpretation of the interaction. In this Appendix we give 
a complimentary physical interpretation of the method to make sense of our momentum 
dependent interaction described in Paragraph 2. For simplicity we restrict ourselves to the 
2-particle Hamiltonian h in equation (2). 

In the main text we gave a formal argument converting the interaction in the Hamil- 
tonian h to the boundary conditions in equation (3). It is interesting to note that, in 
doing this, we have specified a regularization procedure, i.e., given a consistent prescription 
avoiding divergences which would occur in a naive treatment of the singular interaction. 
Indeed, naively the action of /i on a wave function x{x) is {hx){x) = — x"(a;) +4:XS'{x)x'{0), 
but from our discussion in Paragraph 2 it is clear that h is also defined on wave function 
which are discontinuous at x — and with x'(0) therefore undefined. The above-mentioned 
regularization procedure amounts to replacing the ill-defined derivate at a; = by the well- 
defined average of the left- and right derivatives at x = 0, x'(0) ~^ \x'i^~^) + x'{~^~^)]/'^- 
To see that this eliminates a divergence it is instructive to re-derive the bound state energy 
using Fourier transformation. The Fourier transform of hx — Ex can be written as 



where the r.h.s. comes from the interaction with the factor cos{eq) providing the regular- 




(Al) 
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ization and the hat indicating Fourier transform. Computing from this x{k), multiplying 
with k cos{ek) and integrating we get the following self- consistency relation, 



l = 4Alim/^^fi^ (A2) 
e^oJ^2n q^ + \E\ ^ > 

where we used that the bound state energy is negative, E = —\E\. Obviously, without the 
factor cos(eg') the integral on the r.h.s. is linearly divergent, but with this factor we obtain 
the well-defined result 1 = —2X^/\E\, which for A = — |A| has one solution. It is easy to see 
that this yields the same value for the bound state energy and the same bound state wave 
function which we obtained by a different method in Paragraph 2. 



Appendix B. Non-relativistic limits: Details. In this Appendix we give more details 
about how to derive the non-relativistic limits of the Thirring model (Appendix B.l) and 
0f^i-theory (Appendix B.2) discussed in the main text. 

B.l Thirring model. The Dirac Hamiltonian in Eq. (6) in Fourier space is 

with ^2'* = V'±'*(^) ci^iid the hat indicating Fourier transform. It is diagonalized with the 
following canonical transformation. 



where 



Mk) = a±{k)^+{k) ± a^{k)if4k) (B2) 

a±(fc) = yi(l±^, Ek = Vimc^y + {kef- (B3) 
This yields 

Ho^ Jdk: ([Ek - Eo]^l{k)^+{k) - [Ek + Eo]^l{k)^ -{k)^ : . (B4) 

Expanding this in powers of k/me and transforming back to position space one obtains 
the equations given in (8) jf in the main text. Transforming the interaction in equation 
(7) to Fourier space, inserting the equations in (B2), and ignoring the terms involving the 
negative energy fields ^^"^ we obtain 

^it = ^ J dki---dkiS{ki-k2 + h-k^)v{ki,...,ki) ■.^\ki)^{k2)^\h)^{k4) : (B5) 

with the interaction vertex 

v{ki, . . . ,k4) = ^(a+(/ci)a+(fc2)a-(fc3)a-(^4) + a+{k3)a+{k4)a-{ki)a-{k2) 

— a+(/c3)a+(/c2)a_(/ci)a_(A;4) — a+(A;i)a+(/c4)a_(A;3)a_(A;2)) (B6) 
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which we (anti-) symmetrized using the CAR. Expanding this in powers of 1/mc we obtain 

v{h, ...,h) = j^^ih - h){k2 - k,) + 0((mc)-3). (B7) 

Inserting this into equation (B5) and transforming back to position space we obtain the 
interaction term in the non-relativistic Hamiltonian given in equation (9). 

B.2 (j)\^i-theory. This model can be formally defined by the Hamiltonian Ti^ = T-Cq + 



int 



with the free part TCq given in equation (18) in the main text and the interaction 

ng, ^QB j ■.(t>'- (B8) 



with boson fields and 11 as defined after equation (18) in the main text, m > the mass, 
and qb the coupling; the dots indicate normal ordering to be specified below. The free 
boson Hamiltonian in equation (18) can be diagonalized in the usual manner, 

n(x) = --J - 4t(A)e-""') (B9) 

with Ek as in equation (B3) and the <l>^^'' the Fourier transform of non-relativistic boson fields 
$(t) obeying the CCR [^{x),^\y)] = S{x - y) etc. This yields = J dk E^^^ {k)^k) 
where, at this point, normal ordering is defined with respect to the non-interacting vacuum 
|0) obeying $(a;)|0) = 0. Expanding in powers of 1/mc and transforming to position space 
we get 

Ho ^ J dx<^\x)[mc^ -dl/2m + 0{{mc)-^)]<^>{x). (BIO) 

To lowest non-trivial order in 1/mc the first equation in (B9) reduces to equation (22). In- 
serting this into the interaction in equation (B8) we get five terms, but only one of them com- 
mutes with the particle number operator N = J dx $'l"(x)$(x), namely Ggs/i'^.m)'^ J dx : 
[$'''(a;)$(a;)]^ :. The other terms describe processes where the particle number is changed, 
and since the creation of particles requires an energy larger than mc^ (according to equation 
(BIO)) all these processes can be ignored in the non-relativistic limit where mc becomes 
large. ^ Thus the non-relativistic limit of -theory can be described the Hamiltonian 

KVrei = jdx ^^\x){-dmx) + ^ : [^HxMx)r : (Bll) 

which, for 2m = 1 and 3(?B/2m^ = cb, is the second quantization of the ID boson gas 
Hamiltonian given in equation (13). 



^While this is physically plausible, we do not know a convincing mathematical argument to justify this 
simplification. We therefore regard this step as the weak link in our chain of arguments relating the SG 
model to the ID boson gas. 
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